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Abstract 

The AS" = 2 transition b ssd is extremely small in the standard model, proceeding only 
via box diagrams. Such decays are thus an appropriate ground for searching new physics. We 
present a study of several two-body and three-body decays of the B^. meson belonging to this 
class. Our calculation covers the minimal supersymmetric model with and without 7^-parity and 
a two-Higgs-doublet model. The current limits on the parameters of these models allow for rather 
large branching ratios for several two and three-body hadronic decay modes of Be, up to the 10~^ 
range in the minimal supersymmetric model with 7^-parity violation. 
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I. INTRODUCTION 



Rare decays of B-mesons are considered to be one of the promising areas for the discovery 
of new physics beyond the Standard Model(SM) [l|,Q,Q. This is based on the expectation 
that virtual new particles will affect these decays, in particular in processes induced by flavor 
changing neutral currents; such processes are suppressed in the SM since they proceed via 
loop diagrams. This venue, typified by transitions like b s{d)'~f, b s{d)ll has been 
investigated intensively in recent years 

urn. 

A related, though alternative approach is the search for rare b decays which have ex- 
tremely small rates in the SM, and their mere detection would be a sign for new physics. 
Several years ago, Huitu, Lu, Singer and Zhang have suggested the study of ^ ssd 
and b — > dds as prototypes of the alternative method. This proposal is based P] on the 
fact that these transitions are indeed exceedingly small in the SM, where they occur via 
box diagrams with up-type quarks and W^s along the box. The matrix elements of these 
transitions are approximately proportional in SM to A'' and (A is the sine of the Cabibbo 
angle), resulting in branching ratios of approximately 10~^^ and 10~^^ respectively, probably 
too small for detection even at LHC. Further discussions on these (AS* = 2) and (AS' = — 1) 

n n 

Huitu et al. have then investigated |3, 1^ the b — > ssd transition in several models 
of physics beyond SM, namely the minimal supersymmetric standard model (MSSM), the 
minimal supersymmetric standard model with 7?.-parity violating couplings (RPV) and Two 
Higgs Doublet models. Within a certain range of allowed parameters, MSSM predicts a 
branching ratio as high as 10"'' for b ssd. On the other hand, Higgs models may lead 
to a branching ratio at the 10~^ level. Most interestingly, in RPV the process can proceed 

n 

as a tree diagram and the then existing limits ^4] on A-^;, RPV superpotential couplings did 
not constrain at all the b — > ssd transition. 

In Ref. fl th. hadro„,c decay fi- ^ K-R-.-, proceeding either directly or through a 

K*^, has been selected as a convenient signal for the b ssd transition. A rough estimate |^ 
has shown that the semi- inclusive decays i^"i^'" + (nonstrange) may account for about 

1/4 of the b — i> ssd transitions. The OPAL collaboration has undertaken the search for this 



decay establishing 



10| the first upper limit for it, }3n{B- K-R-n^) < 1.3 x 10"^ This 



result then imposes more severe limits on the combination of A^ couplings which appears 
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V obtained from other processes. 



12| the above hmit to BTZ{B 



in the amphtude for h ssd in RPV, than previousi 
Very recently, the BELLE experiment has reduced IllL 
K^K^n^) < 2.4 X 10~^. In an experiment planned [13[ by ATLAS at LHC, a two orders of 
magnitude improvement is expected. 

In order to use these results for restricting the size of the b ssd transition one needs 
also an estimate for the long distance (LD) contributions to such (AS* = 2) processes. A 
calculation ^| of such contributions provided by DD and by Dn intermediate states as 
well as those induced by virtual D and vr mesons leads to a branching ratio BTZ^^{B^ — >■ 
K^K^Ti^) = 6 X 10^^^, only slightly larger than the short-distance result of SM for this 
transition. Thus, this decay and similar ones are indeed suitable for the search of new 
physics, the LD contribution not overshadowing the new physics, if it leads to rates of the 
order of 10^^° or larger. A survey 15] of various possible two-body {AS = 2) decays of B 
mesons to VV, VP, PP states has sing led out B- K*-K^ and B' R-R^ as the most 
likely ones for the detection of the presence of RPV transitions at the 10~^ level. 

In the present paper we study the (AS* = 2) transitions of the B^ meson. This double- 
heavy meson, which has been firstly observed Q, 3| by CDF Collaboration a few years ago, 
can decay weakly by b-quark decay with the c-quark as spectator, by c-quark decay with the 
6-quark as spectator or via an annihilation process. There is already a fairly large number 
of theoretical papers describing the various leptonic, semileptonic and hadronic exclusive 
channels to which B^ may decav in the Standard Model. These are well summarized in a 
recent review by Kiselev jl8l |. 

At the forthcoming LHC accelerator one expects about 5 x 10^'^ Be events/year, at a 
luminosity of 10^^ cm~^s~^ ^]. Hence, the search for rare decays of B^. exhibiting possible 
features of physics beyond SM will also become possible. The effect of such physics on 
radiative decays of i?c, as caused for example by c-quark decay via the c — m + 7 transi- 
tion has already been investigated HQ- Here, we address the effects of new physics on 
rare 6-decays caused by the (AS* = 2) transitions, which we already mentioned to be very 
rare in SM Specifically, we calculate two-body decay modes B^ R , 

D*-R^,D-R*^ and D-R^ as well as three body modes B' D-R-n+, D*-R~t^+, 
D-D*-D+, D-D;D*+, D-D-D+, D^R^R- and D*^R^R~. We expect these modes to 
be most likely candidates for the experimental observation. Our calculation shows that 
while the SM prediction for the branching ratio for these exclusive channels is 10"^'^ or less, 
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some of the {AS = 2) decays can have branching ratios as high as 10 ^ — 10 ^ in MSSM 
with RPV. 

In the second section we describe the framework used in our calculation and the beyond 
SM models which we consider. In Section 3 we give the explicit expressions for the transition 
matrix elements. In Section 4 we discuss the results and the prospects for checking the 
"beyond SM physics" in Be decays. 



II. FRAMEWORK 



The Standard Model analysis describing the b — > ssd transition leads to the effective 
Hamiltonian: 

= CsM (sLYdL) [sLl^hL] , (l) 

where we have denoted 



C: 



1 G 



SM 



2 27r2 



VuVU 



m 



w 



m 



w 



1 ' 2 

ml m^ 



(2) 



with functions f{x) and g{x,y) given in ^. Taking numerical values from PDG 2^ and 
neglecting the CKM phases, one estimates \Csm\ — 4 x 10"^^ GeV~^. 

Next we turn to several models beyond SM; the minimal supersymmetric SM, the minimal 
super symmetric SM with 7^-parity violation and a two-Higgs-doublet model. The largest 



contribution found in the minimal supersymmetric extension of the standard model 
for the corresponding factor 



gives 



C 



MSSM 



216m| 

a 



2AxMx) + mfeix) 



with as{rnw) = 0.12, x = m-fm"-, and the functions fei^x) and fei^x) given 



in 0- 



(3) 



The 



couplings Sfj parametrize the mixing between the down-type left-handed squarks. The value 
of 6i2 is determined from the — K^ mixing {2^ and is currently bounded by —mxg = 
3.49 X 10"^^ GeV {2^. We present our calculation choosing x = 1 and 8; using results 
from [2^ we estimate the absolute value of 612 to be below 3(9) x 10^^ at average squark 
mass mj = 350(500) GeV and x = 1(8). The strongest bounds on ^23 come from the 
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radiative b ^ s'y decay [23|, |25[. These studies give at x = 1 and for = 350 GeV 
the stronger bound on |523| ^ 0.4 which results in \Cmssm\ ~ 5 x 10^^^ GeV^^. Taking the 



same values as in 
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x = 8, mj = 500 GeV, leads to \Cmssm\ ^ 3 x 10"^° GeV"^ These 
updated values for Cmssm are somewhat smaller than those used in P. IfS^. 

When including i?-parity violating interactions in the MSSM, the part of the superpoten- 
tial which is relevant here is ly = X'^jf^LiQjdk, where i, j and k are indices for the families and 
the L, Q and d are superfields for the lepton doublet, the quark doublet and the down-type 
quark singlet, respectively. Using the notation of j3| and j^, the tree level Hamiltonian has 
the form 

3 



= -^^^[K32K2lisRbL)isLdR) 
n=l ^'1 

+ Kl2K23isRdLKsM]. (4) 



The QCD corrections were found to be important for this transition j27|. For our purpose 
it suffices to follow retaining the leading order QCD result /qcd — 2, for mn = 100 GeV. 

The most recent upper bound on the specific combination of cou pling s entering (j3]) can 
be obtained from BELLE's search for the K^K^tx^ decay (ill. Il2j|. The result based 

on the explicit calculation of the amplitude B~ — > K'K^tt^ j28| 



aRPV _ IqCD rw 

n=l 



and experimental upper bound from which corresponds to branching ratio BTZ{B 



K-K-TI+) < 2.4 X 10^6 is yEn=i \Kz2K2i\ + \K2i^2L< 3.7 X 10-^ Note that this 
is a somewhat larger value than previous estimates [l^ |2^ based on the assumption of 
V{B- ^ K-K-TT+) ~ l/4r(6 ^ ssd) J. 

The decay b —>■ ssd has been investigated using two Higgs doublet models as well . 
These authors found that the charged Higgs box contribution in MSSM is negligible. In two 
models with charged Higgs branching ratios are very small ^ for reasonable values of cot (3. 
Recently it has been also verified that the chargino contribution in MSSM to this process 
is indeed smaller by an order of magnitude than contributions calculated in (41. On the 

n 

other hand, THDM involving several neutral Higgs bosons [22|| could in principle give larger 
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inclusive decay rates for b ssd. The part of the effective Hamiltonian relevant in our case 
is the tree diagram exchanging the neutral Higgs bosons h (scalar) and A (pseudoscalar) 



TH 



2 ^sb^s 



-^{sd){sb) V(s75(i)(s756) 



m 



(6) 

with the coupling defined in j29|] as a Yukawa coupling of the flavor changing neutral 
current transitions di ^ dj. In our estimation we use the bound \^sbC,sd\/nT''H > lO^^*' GeV~^, 
H = h, A which was obtained in by using the Auik limit on ^bd/'m'H and assuming 

lU/mnl > 10-^ 

In order to estimate the amplitudes for the various AS* = 2 exclusive decays of B^, one has 
to calculate the matrix elements of the operators appearing in the effective Hamiltonian. We 



rely on the factorization approximation 



30 



31 



32j | which requires the knowledge of matrix 



elements of the current operators or density operators. More sophisticated methods such 
as perturbative QCD calculation have been employed in the 5° case However, 

in the case of nonleptonic weak decays of Be mesons such approaches have not yet been 
developed. We use the standard form factor representation of hadronic current matrix 
elements 



31 



m 



35 



. For the matrix element between two pseudoscalar mesons 



we employ 



{P2{P2mYqr\PliPl)) 



m 



m 



(7) 



where = {pi — P2)^ and the two different parameterizations are connected via -^0(9^) = 
f+{q^) + f-{q^)q^/{'fnpi — '^P2) and Fi{q^) = f^{q^). The matrix element between a pseu- 
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doscalar and vector meson we decompose as customary 



mp + my 
—i{mp + my 



e ■ q 



2mv 

Mq') 



(mp + my) 



iPl+P2) 



mn — m 



A2{q' 



-Fy{q')e^^^^e*y,p,^P2p + iF^{q^)e'^ 
+iF^{q^)e*y ■ p,{pi + p^r + ^F\*y ■ p,q>^, 



(8) 



and again the two different parameterizations are connected by Vlq"^) = Fy{q^){mp + my). 
Ao{q') = {l/2my){-F^{q') + F^{q'){ml-ml) + F^{q^)q'), A,{q') = F^{q') / {mp + my) 
and A2{q^) = F^{q^){mp + my). We also use the following decay constants 



(P(p)|g,7Vg.|0) = ifpp', 
{S{p)\q,Yqm = fsP", 
{y{(^v,p)\qjltiqi\0) = e*^y(p^), 

Finally, for the calculations of the density operators we use relations 

(9°(s7q,6) = i{mb — ms)sb 
d'^is-falbb) = i{mb + ms)sj5b. 



(9a) 
(9b) 
(9c) 



(10a) 
(10b) 



In our calculations we need the oi*^ transition form factors /±, Fy, F^ and F^. 

Since heavy quark effective theory (HQET) is not directly applicable to the decays of the 
Be meson, we assume pole dominance for these form factors |3ll . |36| : 

m 



F{q' 



q^ /m^ 



(11) 



pole/ 



and take numerical values for -F(O) and mpoic from from QCD sum rules calculations 
(see Tables HI and ITT]) . For the F^'^ and F^'^, which are needed in the three body decay 
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TABLE I: Numerical values of Be — > Dg transition form factors at g'^ = taken from 



36| 



/+ /_ Fy(GeV)-^ F^^iGeY) {GeVy' (GeV)- 



0.45 -0.43 0.24 


4.7 -0.077 


0.13 


TABLE IL Pole masses used in B^, — 


^ Dg transition form factors taken from j 


3- 


UiGeV) /-(GeV) Fy(GeV) 


Fo^(GeV) F4*(GeV) 


F^(GeV) 


5.0 5.0 6.2 


oo 6.2 


6.2 



modes involving K and 7r+ mesons we use results from js?! 



?K7rf2^ _ '^9VKi892)9K 
1 W 



T - mj^(892) + WT^K{892) 
2fi'yx(892)fl'ii-'(l - Q'^/"^/C(892)) 



- ^K(892) + W ?^ri^(892) 



"1^(892) 

_l (f f K (1430) 9 SK (1430) q^n 

(ml - ml) q2 _ m^^.^g^) + 2 V^r^(i43o) 

In our numerical calculations we use the following values: Jk = 0.162 GeV, /a'o(1430}_^ 
0.05 GeV, qk* = 0.196 GeV^, 5'v'x(892) = 4.59 and 1430) = 3.67±0.3 GeV taken from |37l |. 

In the three body decay modes involving pairs of D and Dg mesons, we also need the 
form factors for the D^*^ D^*'> transitions. These are not available in the literature and 
we calculate them by utilizing HQET, including the light scalar meson interactions with 
heavy mesons as presuming the main contributions from 



exchange of light meson resonances {K^, K*{892) and -ftr*°(1430 



mesons are described by the HQET chiral Lagrangian 



m. 



Interactions of heavy 



+ipTr [Hv^iVf" - p^)^ - ^Tr [Hall] + ... (14) 

We use standard heavy meson field notation H = 1/2 + 1) {P^'y^ — -P75), where P and P* 
create heavy pseudoscalar and vector mesons respectively, and exponential representation of 
light pseudoscalar meson fields ^ = exp{iA4/ f), where JH is the light pseudoscalar meson 
chiral SU{3) nonet matrix, so that = 1/2(^^8,,^ + ^8,,^^), = l/2({ta^^ _^5^^t). Light 
vector fields are represented as = igvPij./V2, where Pfj, is the light vector meson chiral 



SU{3) octet matrix. Finally light scalar mesons are introduced through the a = a/2/3(J 
field, where a is the light scalar meson chiral SU{3) nonet matrix. The ellipses indicate 
further terms involving only light meson fields, chiral and I/itih corrections. The light 



vector meson coupling {(3) was found to be consistent with zero 42]. Consequently we 



neglect the contribution of the vector fC*(892) resonance in our calculations. In the HQET 
approximation in the leading order in heavy quark mass and chiral perturbation /_,_ " , Fy " , 
F^'^\ F^^''^, F^^'^*, F^^*^ and F^^'^* are all found to vanish, so the only contributions 
come to the F:"^^*^ and form factors 

fDsD^^2^ ^ {9j^)fK{um^mDmD ^^^^ 

'A'(1430) 



(16) 

In the numerical calculation we use g^, =~ 3.73 and qa — 0.59 0]. The rest of parameters 
are taken from PDG 

III. AMPLITUDES 
A. Two-body decays 

First we investigate two-body decays of B~ mesons which a priori are the promissing 
candidates for the experimental studies. 

1. B- ^ D*~K*^ decay 

For the analysis of pseudoscalar meson decay into two vector mesons it is convenient to 
use partial helicity formalism (see, e.g., The amplitude can be decomposed in the 

product of the effective Wilson coefficient and the matrix element of the operator. Within 
SM we have 

M{B; ^ Dl-R*"^) = ^Csm{D:-K*''\0\B-), (17) 

where O = (s7''(l — 75)rf)(s7''(l — 75)6). Using factorization one finds {D^^ K*'^\0\B^) = 
— {D*^\{s'~f^b)\B^) X (A'*^ I (57^756?) 1 0) which by using helicity formalism leads to the follow- 



9 



ing expressions for amplitudes: 

A^oo 

M± 



-;CsMgK'{mB, + niD*) 



(18) 



x[aA^{ml.)TlV{ml.)l (19) 
where a = (1 — — t^)/2rt, /3 = k"^ /{2rt{l + r^)), 7 = A;/(l + r^) with r = mD^/mK", 



t = rriK* /friB^ and A; = ^/l + r'^ + — 2r'^ — 2t'^ — 2rH'^. The decay with is then 

|A^oor + |A<++r + |A<-|']. (20) 



r(i?e 



b1 



Within the MSSM the result in (^01) can be reused by replacing Csm with Cmssm- On 
the other hand RPV and THDM models cannot be tested in this mode when factorization 
approach is used due to the vanishing {K*'^{pk)\{sj^{1 — 75)6?) |0) matrix element. 

2. B- D*-K^ decay 

In SM and MSSM factorization gives for the O matrix element {Dl-K^\0\B-) = 
{K^\{sr{l - l5)dm{D*-\{sr{l - I^MB;) = -2tmo*fKAo{ml)e}, ■ pk- The decay 
width becomes 



\p\ 



SM 



j2\{D:-K'mB;)\' 

pol. 

where J2poi sums over all polarizations of D*~ meson. One finds J2poi \^*d ' PkI"^ = 
A(m| ,m|^,m|,,), where X{a,b,c) = {a? + 6^ + — 2{ah + he + ac))/Ac. RPV can also 



(21) 



occur in this decay. The decay amplitude is 

M{B- ^ D'-K") 



E 



^'^^^^'^'' {d:-k^\Or,\b;) 



4 



(22) 
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where we have denoted Or^ = — 75)&) + 75)^) and Or^ = (s(l — 75)6?) (s(l + 75)&). 
Factorization yields 



{D:-K'\OnAB:) = {D:-K'\Or,\B-) 

(^°|(.-(l-75M)|0)(Dr|(s(l + 75)&)|S.-) 

(2mz3*e;^ -Pk) ^o("^|-), 



m^fK 



{nis + rrid) (m^ + rrib) 



(23) 



giving for the decay width 



\P\ fqCD 
7rml mi 



/* l'^ I I \' \'* 1^ 



pol. 



Similarly, the THDM gives for the amplitude of this decay 



M{B- ^DrK') = --UU 



X 



^ -{D:-K'\OT^B-) - 1^{d:-k'\ot,\b-) 



nit 



(24) 



(25) 



where C^i = {sd){sb) and = (s75C?)(s756). Only the second amplitude turns out to 
be non- vanishing and gives the same result as ()23|) . Because of the specific combination 
of products of scalar (pseudoscalar) densities, this is the only two body decay which has 
non-vanishing THDM amplitude within the factorization assumption. 



3. B- D~K^* decay 

For this decay mode the matrix element of the operator O is determined to be 
{D^K*^\0\B-) = {K*^\isr{l - l5)d) \0){D-\{sril-l5mB-) = 2gKF,{ml.)e*^ ■ p^. 
On the other hand both RPV and THDM amplitudes for this mode vanish, due to the 
vanishing of the matrix element of the density operator for K*^ state. 
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4- B- ^ DJK^ decay 



The matrix element for the operator O becomes in this case {DgK^\0\B^) — 
{K^\{sY{l - i,)dm{Dj\{sY{l - l^)h)\B-) = ifKFo{ml){ml^ - mlj, while for the op- 
erators Or^ and Or^ we get 

{D:K'\On,\B:) = -{D:K'\Or,\B-) 
= (XO|(.-(l-75)rf)|0)(i?;|(«-(l + 75)&)|S;) 

(26) 

THDM does not contribute to this mode. 
B. Three-body decays 

Next we study three body decay modes of B~ which are also possibly good candidates 
for experimental searches. 

1. B~ — ^ DjK~'iT'^ decay 

The required matrix element in SM and MSSM is here 

{d:k-7:^\o\b:) 

= (X-7r+|(s7'^(l - ^,)dm{D:\{sr{l - ^,)b)\B-) 

= /r(^)/r(^)(^ - ^) + ^(^)(^?. - o 

- ml) + ^(5)/5'^(5)(-5), (27) 

where s = mig = = (ps - Pd)^, u ^ ml^ ^ {pb - p^)^ and t = mfg = {ps - Pk)^. We 
obtain the decay rate by integrating the transition matrix element over the whole Dalitz 
plot 



X / dm^2 / dml^ 



(27r)3 32m|^ 

1 ' ' 

-CsM{D:K-7r+\Oi\B; 



(28) 
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where {ml2)max = {jriB^-m^Y, {rnl^)^^ = {niD + mKY, while 
{ml^)max = {E; + Elf - {^Ef -ml - ^ Ef - mlf and {ml^)min 
{E* + Elf - {^Ef -ml + ^Ef - mlf. E* = {m\^ - ml + ml)/2mu and 
El = {ml^ — ml — ml) /2mi2. We also get nonvanishing contribution in the RPV 
model from the matrix elements 

{D^K-'n+Pn^B-) = -{D; K-n^\On,\B;) 

= {K~rc+\{s{l-^,)d)\0){D:\{s{l + ^,)b)\B;) 

- F^^(c)F^-(c) ^""^^ ~ - (291 

" ° ° ^ ^ (-m. + m,)K-m.) " ^^^^ 

Finally, THDM contributes to this decay mode with the matrix element for the Ot2 operator 
the same as 



2. B- D*-K~TT+ decay 

In SM and MSSM the matrix element of the operator O can be factorized as 

{D:-K-n+\0\B-) = {K-n+\{-sY{l-l5)dm 

x{D:\{sr{l-l,)b)\B;). (30) 

The actual expression of the amplitude A in terms of the corresponding form factors turns 
out to be rather long and complicated. This is is even more so for the squared amplitude 
Spoi l-^P' once we have summed over polarizations of the D*~ meson. Thus, we only give 
for this mode in SM and MSSM the final numerical values of the branching ratios in the 
results section. While the THDM does not contribute to this decay mode, the RPV model 
also gives a non-vanishing contribution from the matrix elements 

{D:-K-n+\OnAB;) = {D*- K-n+\On,\B;) 
= {K-rr+\{s{l - ^,)d)\0){D:\{s{l+j,)b)\B;) 



-2i 



<^(s)mB| - Af^(s)(mB, + mo*) 



Pb ■ en 



><F^^{s)- ^^-^ -, (31) 

[-ms + md){mb + ms) 

while summing over polarizations of the D*~ meson gives Ylpoi \Pb ' ^-Op = -^("^5^5 "5, ''^i)*)- 
Note that expression (|HT|) involves a slightly different combination of B^ —>■ D*~ form 
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factors than for example the analogous expression ()23|) . This is because we can assume 
2mDj^o("^l:) — i^Bc + ^D*)Ai{mj^) — {rriB^ — rnD*)A2{m'j^) at small kaon masses, while 
the same is not true for all values of s in the kinematic range. 

3. B- D-D*-D+ decay 

Here we get two contributions to the decay amplitude. The matrix element for the 
operator O becomes 



The second term turns out to be relatively suppressed due to small values of /i^*(i43o) 
and (yfTT, in the F^"^ form factor while averaging over polarizations in the first term gives 
XI \^*Ds ■ ~ -^(Q'^) "^i)) "^dJ- Similarly in the RPV model we only keep the first of the two 



possible amplitudes {Dj D*- D+\Or,\B;) = {D;D*-D+\Or,\B-) ^ {D;\s{l - -f^)b\B-) x 
{D*-D+\s{l + 75)rf|0), so that 



{D;d:-d^\0\b;) 

= {D;\-sY{1 -l')b\B;){D:-D+\-s^,{l -^,)d\0) 
+ {D:-\-sY{l-l')b\B;){D;D^\-s^,{l-^,)d\0) 

+ {ml-ml)F^^^\q')F,^^^{q')e^,- -q 



(32) 



/ 2 2 \ 2 



(33) 



4- B- 



D-D-D*+ decay 



The matrix element of the O operator becomes in this case 



{D;D-D*+\0\B;) = {D;ipnJ\rsYil-l')b\B:{pB)) 
x{D;ipnjD*^ipn,)\rs^,{l-l,)d\0) 



(34) 
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so that Epoi. = \Csm/M' X f{t){2f{s){ml. - mlj + f{t) (t-u)), where f{x) = (m|^ 
— mlj X F^^'^ [x)Fq^[x). In the RPV model we get for this mode equal matrix elements 
of operators Or^ and Or^ 

= {D;{pj,J\s{i-^')b\B;{pB)) 

x{D;{pD,)D*+{pD,)m + -f5)d\0) 

+ {Pdi^Pd2), (35) 

so that Epoi. \M\' = ifQcn/mlfxij:^ \K32K21? + \K2iK32\') x 9im9{s){ml, -m^J + 
g{t){t — u)), where g{x) = {m%^ — m|,J x F^^^^ {x)Fq^{x) x x/ (m^ + md){mb — rus). 



5. Dg Dg decay 

The matrix element of the operator O is calculated to be 



{D:D:D+\0\B;) = {D:{poJ\sr{l-j')b\B;{pB)) X {D:{pn.)D+{pn,)\sj,{l-j,)d\0) 

+{s ^ t) (36) 
RPV model also allows for this mode 



{D;D;D+\OnAB;) = {D;D;D+\On,\B;) 

= {D;ipn,)\s{l-^')b\B;ipB)) 

x{D;{pD,)D+{pD,){q)m + ^5)d\0) 

+ {PDi ^PDi) 



X- 



Finally, the THDM contributes to the amplitude of this decay through the operator Oti 
giving the same result as (j^Zj). 
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6. B- K^D^K- decay 

Also this transition has nonzero amphtudes in all models considered. We analyzed this 
decay amplitude with the assumption of D*~ vector meson resonance exchange and used 
HQET for the evaluation of the D*~D^K vertex. The matrix element of the operator O 
turns out to be 

{K'^D°K-\0\B-) 

= {D:-iq,en\sY{^-l')b\B;iPB)) 
x{K\pK,)\rs^,{l-^,)d\0) 

xtDj^^iq) X {D\pD)K-ipK2)\Cs\D:~{q,en), 

(38) 

where D^^^ (q) = {—g^u + QpL(lu/^'D*)/{Q^'-^'D* + i^{Dl)mj:,*) is the D*~ meson propagator 
while {D^{pD)K~{pK2)\C-s\Dl~'{q)) = {2gA/ fx) x y/mD^niD x ■pK2 is the strong transition 
matrix element calculated in HQET. The amplitude thus becomes 

MiB; ^ K'D'K-) = -CsM2ifKAl''{m\) 

X 7i — ; 7 r 7j—t2-—^/mn*mn. 



(39) 



Similarly in the RPV model we get 

= {k'^d''k-\Ob.,\b;) 

= {D:-{q,en\sil-^')b\B-{pB)) 
x{K\pKi)m+l5)d\0) 

xiD'^Jiq) X {D\pn)K-{pK2)\Cs\D:{q,n) 

2ifKmlA^^{ml) 
(mfo + ms){md + m^) 

{s + ml,){u-2mj^) 

2mjj,{s + tr{D*)mD* -mf),) f^^ 



(40) 



Finally in the THDM only the operator Oti has non-vanishing amplitude and gives the same 
result as (gn)). 
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7. B- K^D*^K- decay 



This decay model was analyzed by presuming pseudoscalar meson resonance ex- 
change. The matrix element of the operator O was found to be 

{K^D*'^K-\0\B-) 

= {D;{q)\sY{l-^')b\B;iPB)) 
x{K%pKi)\s^,{l-l,)d\0) 

x^A^^(g) X {D*\pn)K-ipK2)\Cs\DM), 

(41) 



where we denoted the pseudoscalar meson propagator with A^°(g) = l/(g^ — mj-ij. 
Evaluating the strong interaction vertex with HQET we get for the amplitude 



m 



g2 - mf,^ 



,9 a 



xi2—y/mD,mD*e*j^ ■ Pk2- 
Ik 



(42) 



RPV model also contributes to this decay mode, giving for the Ori and Oji2 niatrix elements 

{k^d*^k-\Or,\b;) = -{k^d*^k-\Or,\b;) 



{D;iqMl-^')b\B;ipB)) 
x{K°{pKi)m+l5)d\0) 
x^A^=(g) X {D*%pj,)K~{pK2)\Cs\D;{q)) 



'K) 



{nib - rns){rnd + m. 



m 



X ■ 



Br 



— m 



-^i2 — ^mDsfnD'e^ ■ Pk2- 



K 



IV. RESULTS AND DISCUSSION 



(43) 



The usefulness of AS* = 2 decays o 



discussed in several publications 



' B mesons in the search for new physics has been 
IflaBQQ.FVo,„*e.„de.c„„.de..eds„ 
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far it appears that these decays are particularly relevant in the search for supersymmetry, 
with and without 7?.-parity violation. In the present paper we have undertaken a similar 
investigation with regard to the AS = 2 decays of the meson. 

The calculation of matrix elements has been performed in the factorization approxima- 
tion. In the 5-meson decays this works in a good number of cases, while in other cases a 
more sophisticated approach is needed (for a recent review see Q]). In this first calcula- 
tion we considered that factorization approximation is sufficient for obtaining the correct 
features of the decays of the various channels considered. An exception is the case in which 
matrix elements vanish as a result of factorization, which in a better approximation can be 
improved. 

The results obtained in the MSSM framework depend on the values of the 6fj parameters 
of the mass-insertion approximation which we use. The constraints on these parameters 
have been improved in recent years Q| vs. the values which were used in the original 
calculation of the b ssd. For the 7^-parity violating MSSM we use the new limits on 
the A'y couplings of the superpotential, which we obtain by the use of the new limits on 
the — » K^K^-K^ decay from BELLE ^] and the explicit calculation of this decay of 



reference 

In Table IIIII we summarize the results of our calculation for the various possible decay 
modes considered. As expected, the SM results for all these AS = 2 channels are extremely 
small. For the models beyond SM, the values presented are the upper limits presently allowed 
by the constraints available from other processes. For MSSM we present the results for the 
two extreme values of x = 1 and 8. We note that MSSM with 7^-parity violation is rather 
poorly restricted at present for some of the decays, in particular the decays B^ — D*^K^, 
B- D*-K'n+, B- D-K~TT+, B' DjD*-D+, B' D;D-D*+, which are in the 
10^^ — 10^^ range. All these decays should be looked for, when sizable samples of B^s will 
be available. 
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Decay SM MSSM (.x = 1) MSSM (x = 8) RPVM THDM 



Be 




B>s ^0 


O X 


in-14 


X 


1 n-13 
iU 


Z X iU 












Be 




Ti*— f?' 


X 


in-l4 


O X 


in-l3 
iU 


1 xy 1 n — 9 
i X iU 


i X 




-5 


10 


-10 


~r> — 

Be 




B>s 


D X 


in-l5 
iU 


4 X 




o xy 1 n — 10 
z X 10 












Be 




7~i — f?' 


4 X 


1 n-15 
lU 


Z X 


1 n— 14 
ID 


n 1 n — 11 
9 X 10 


1 X 


lU 


-6 






Be 




Ug K vr^ 


Z X 


1 n-14 


i X 


1 n-13 
iU 


A -.y in — 10 

4 X 10 


D X 


1 n— 

10 


-5 


10 


-10 


Be 




JJg K TT^ 


D X 


in-14 


4 X 


1 n-13 


1 x^ 1 n — 9 
1 X lU 


n x^ 
9 X 


1 n~ 
10 


-5 






B7 

C 




s s 


3 X 


10-1^ 


2 X 


10-16 


7 X 10-13 


5 X 


10" 


-5 






Be 




D~D-D*+ 


1 X 


10-1^ 


7x 


10-17 


3 X 10-13 


2 X 


10- 


-5 
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D-DJD+ 


9 X 


10-20 


6 X 


10-19 


2 X 10-15 


2 X 


10- 


-7 


10- 


-12 


Be 




DoKoK- 


2 X 


10-14 


1 X 


10-13 


6 X 10-1° 


6 X 


10- 


-6 


10- 


-11 


Be 






2 X 


10-14 


1 X 


10-13 


4 X 10-1° 


5 X 


10- 


-6 







TABLE III: The predicted branching ratios for the AS = 2 two- and three-body decays of B~ 
calculated using the factorization approach within Standard model (the first column), Minimal 
supcrsymmctric standard model for x = I and .x = 8 as explained in the text (the second and third 
column), Minimal supersymmetric standard model extended by TZ- parity breaking (the fourth 
column), and two Higgs doublet model (the fifth column). In columns (2)- (5) the calculated upper 
limits are given as allowed by presently known limits from other processes. 
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